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ABSTRACT 

We analyse a set of galaxy interactions performed by using a self-consistent chemo- 
hydro dynamical model which includes star formation, Supernova feedback and chem- 
ical evolution. In agreement with previous works, we find that tidally-induced low- 
metallicity gas infiows dilute the central oxygen abundance and contribute to the fiat- 
tening of the metallicity gradients. The tidally-induced infiows trigger starbursts which 
increase the impact of SN II feedback injecting new chemical elements and driving 
galactic winds which modulate the metallicity distribution. Although a-enhancement 
in the central regions is detected as a result of the induced starbursts in agreement 
with previous works, our simulations suggest that this parameter can only provide 
a timing of the first pericentre mainly for non-retrograde encounters. In order to re- 
produce wet major mergers at low and high redshifts, we have run simulations with 
respectively 20 and 50 percent of the disc in form of gas. We find that the more gas-rich 
encounters behave similarly to the less rich ones, between the first and second peri- 
centre where low-metallicity gas infiows are triggered. However, the higher strength of 
the infiows triggered in the more gas-rich interactions produces larger metal dilutions 
factors which are afterward modulated by the new chemical production by Supernova. 
We find that the more gas-rich interaction develops violent and clumpy star formation 
triggered by local instabilities all over the disc before the first pericentre, so that if 
these galaxies were observed at these early stages where no important tidally-induced 
infiows have been able to develop yet, they would tend to show an excess of oxygen. 
We find a global mean correlation of both the central abundances and the gradients 
with the strength of the star formation activity. However, the correlations are affected 
by orbital parameters, gas infiows and outfiows, suggesting that it might be difficult to 
determine it from observations. Overall, our findings show that a consistent descrip- 
tion of the gas dynamics and stellar evolution along the interactions is necessary to 
assess their effects on the chemical properties of the interstellar medium. 

Key words: galaxies: formation, galaxies: evolution, galaxies: interactions, galaxies: 
abundances. 



1 INTRODUCTION 

Galaxy interactions are considered a key process in 
galaxy formation since these violent events could be 
able to redistri bute mass and angular momentum ver^ ' 



able to redistri bute mass and angular momentum very 
efficiently (e.g. iBarnes k Hernouistl Il996l : iTisseral l2000l : 
iDi Matteo et al.1 l2008h . As a consequence, chemical ele- 
ments can be stirred up affecting their distribution and 
metallicity patterns such as gradients. There are numer- 
ous obser vational works dedicated to the study of galaxy 
pairs fe.g. ISersic Pastorizalll967l: iTinslev Larson| | 19781 : 
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iBarton et al.ll2000l : iLambas et al.H2003l : iNikolic et alJl2004 



Alonso et all l2006l: IWoods et all l2006l : lEllison et aP l2008l . 
20101 : IPatton et aLllioTlh . But only in the last few years, 
it has been possible to focus on the effects of galaxy 
interactions on chemical abundances. Interacting galax- 
ies, galaxy pairs, ultra- luminous infrared galaxies, and 
merger remnants have been found to exhibit low gas- 
phase metallicity at a given stellar m ass compared with 



RuDke et al. 20081: 


Ellison et al. 200^: Michel-Dansac et al.1 


2008; "Peeples et al. 
Kewlev et al.. (.201( 


|2009l:ISol Alonso et al. 2010). Recentlv. 
J) analysed the metallicity gradients of 
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galaxy pairs in detail confirming the central dilution and 
gradient flattening previously reported by several authors. 
These trends have been interpreted as a consequence of 
low-metallicity inflows triggered by galaxy interactions. In- 
deed, these interacting galaxies are in majority spirals 
which, in the nea rby Universe, hold a gas-phase metallic- 
ity gradient (e.^ IZaritskv et al.1 Il994l : jputil k Rovl Il999l : 
iMoustakas et al.l I2OIO ') . The presence of metallicity gradi- 
ents could explain the decrease of the central abundances 
when low-metallicity inflows are triggered during interac- 
tions. There are also some galaxies in pairs with abun- 
dances departing signiflcantly from the m ass- metallicity re - 
lation such as t hose re port ed by Michel- D ans ae et al.l (|2008h : 
IPeeples et all (|2009l ): fSoTAlonso et al. (201 d). Hence, there 
are still effects of galaxy interactions which demand further 
understanding. 

B y using cos n iologi cal chemo-hydrodynamical simula- 
tions, (Perez et"aI1 (|2006[ ) showed for the first time that the 
interacting galaxies tend to have lower central metallicity 
when compared to isolated galaxies of similar stellar mass. 
These authors showed that the decrease of central abun- 
dances could be explained in terms of gas inflows triggered 
during the interactions. More recently other numerical works 
have conflrmed t hese flndings by using different simulation s 
and hvDothe sis (iRupke et al.ll2010al : iMontuori et aLllioToh . 
Particularly, Rupke et al.l (|2010ah studied the evolution of 
metallicity proflles by using hydrodynamical simulations to 
recreate the major interaction of galaxies with a gas mass 
fraction similar to that used in our less gas-rich simulations. 
However, these authors did not consider new star formation 
activity and the subsequent chemical enrichment of the in- 
terstellar medium as we do in our work. 

In this paper, we use a chemical model grafted onto 
GADGET-2, which allows the description of the chemi- 
cal enrichment of baryons as the merger evolves. Our 
code also includes a self -consistent model of SN feedback 
(IScannapieco et aL I I2OO6I ). so mass-loaded galactic outflows 
can be triggered without including mass-d epend parame- 
ters. T his code includes the chemical model of'Moscon i et al.l 
which describes the enrichment of baryons by Type 
la and Type II SNe. Hence, the simulations analysed in this 
work follow the dynamical evolution together with the chem- 
ical evolution of baryons as the interaction takes place. 

This paper is organized as follows. In Section 2, we de- 
scribe the code and the initial conditions designed for this 
work. In Section 3, we explore wet interactions at low red- 
shift while Section 4 analyses a more gas-rich encounter, in 
order to reproduce typical wet interactions at high redshift. 
Finally, we summarize our flndings in Section 5. 



2 NUMERICAL SIMULATIONS 

We performed a set of hydrodynamical simulations of in- 
teracting pre-prepared disc galaxies of comparable masses. 
The simulation s were run by using a version of GApGET-2 
(|Springell l2005h modifled by 'Scannapiec o et al.1 (EooBI. I2OO6I I 
to include the chemical model of Mosconi et al ] (l200lh and 



a new SN feedback model developed to provide an improved 
description of the interstellar medium (ISM) and the injec- 
tion of energy and metals generated by S Ne. The multiphase 
and SN feedback model developed by IScannapieco et ahl 



(|2005l , l2006h allow the coexistence of gas clouds with dif- 
ferent thermodynamical properties and are able to describe 
the injection of energy from young stars into the ISM pro- 
ducing the self-regulation of the star formation activity as 
well as triggering galactic outflows without the need to in- 
troduce mass-dependent parameters or to change discontin- 
uously particle momentum to start a wind or temporary sup- 
pression of radiative cooling. This scheme has been success- 
ful at reproducing the star formation activity during quies- 
cent and starburst phases and is able to trigger mass-loaded 
winds with strength dependent on the potenti al well of the 
systems (e.g. [Scannapieco et al. 2006; De Ros si et al.ll2010l : 
ISawala et al.ll2011l : Scannapi eco et al.ll2011I V 

Chemical elements are synthesized in the stellar interi- 
ors, and ejected into the ISM by SNe. The chemical model 
describes the enrichment by SNII and SNIa, providing in- 
formation on certain individual chemical elements. We also 
assume that each SN injects an amount of energy into the 
ISM, ^sN (cf Table [U). SNII and SNIa are originated from 
different progenitors with different rates, yields and time- 
scales. We assume that SNII are generated by stars more 
massive than 8 Mp after ^ 10^ y rs. We adopt the yields 
given bv IWooslev Weaver! (|l995l V For SNIa the chemical 
model assumes lifetimes randomly distri buted in the range 
[10^, 1 0^] yrs and yields from W7 model of lThielemann et al.1 
(119931). The SN feedback model is able to generate galac- 
tic outflows powerful enough to transport enr iched material 
outside the galaxies ([Scannapieco et aT]|2008l V Thus, these 
simulations provide us with the evolution of the dynamical 
and the chemical properties of the baryons as galaxy inter- 
actions proceed. 

We run three major merger simulations (1:1) varying 
the orbital parameters (SimI, Simll, and SimV) with the 
same star formation and SN feedback model. We also per- 
formed a run (SimlV) with the same initial conditions as 
Simll but decreasing the energy release by each SN event 
(E^at) and a run (Simlll) with the same orbital parameters 
but without star formation and feedback activ ity. This last 
experiment is comparable to those analysed bv lRupke et al.l 
f2010a). In order to represent wet mergers at low redshift, 
these experiments (from SimI to SimV) were run with a 
20 percent of the disc in form of gas. Finally, SimVI has 
the same orbital, star formation and feedback parameters 
as Simll but the gas fraction in the disc has been increased 
to up 50 per cent in order to simulate a wet encounter at high 
redshift. The parameters of the simulations are summarized 
in Table [T] 

Initial conditions correspond to disc galaxies composed 
of a DM halo, a bulge component (modeled by a NFW and 
a Hernquist proflles, respectively) and an exponential disc, 
with a total baryonic mass of Mb - 5 X 1O^°M0. For the less 
gas-rich simulations, we use 200 000 dark matter particles, 
100 000 stars distributed in the stellar disc and the bulge, 
and 100 000 particles to represent the disc gas component 
with an initial gas resolution of ~ 3 x 10^ M0. For the more 
gas-rich experiment, we re-calculate the number of particles 
so as to keep similar mass resolution. Total stellar masses 
change according to the gas content in the galaxy discs, be- 
ing M* 4.4 X 10^° M0 for the less gas- rich experiment and 
M* - 3.2 X 1O^°M0 for the more gas-rich simulation. A 
gravitational softening of ec = 0.16 kpc was adopted for the 
gas particles, es = 0.2 kpc for the stars, and cdm = 0.32 kpc 
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Table 1. Characteristics of the numerical experiments, i: angle 
between the two total angular momentum of the galaxies. E^^v: 
amount of SN energy release by each event, /gas: gas fraction of 
the discs. 



Simulation 



i n EsN (10^1 erg/s) /g. 



SimI 


45 


0.5 


0.20 


Simll 





0.5 


0.20 


Simlll (no SF) 





0. 


0.20 


SimlV 





0.04 


0.20 


SimV 


180 


0.5 


0.20 


Sim VI 





0.5 


0.50 



for the dark matter. Galaxies live in DM halos with circular 
velocities of 160km/s. All these experim ents are set to have 
the sa me elliptical orbits, according to iToomre Toomr3 
(|l972l l, with pericentre distance of ^ 20 kpc and circularity 
parameter e ^ 0.20. These orbital values are consistent with 
thos e obtained from the analy sis of cosmological simulations 
fe.g. lKhochfar &: Burkert|[2006 V We have kept them fixed so 
that changes in the chemical and dynamical properties can 
only be adscribed to the different angular momentum orien- 
tation. A larger orbital distribution might be of interest to 
analyse possible dependences in the future. 

The gas component is initially assumed to have been 
pre-enriched by the existing stars, following a metallicity 
gradient statistically consistent with observations. We as- 
signed an initial amount of chemical elements to the gas 
particles in order t o reproduce the me an observational gra- 
dient reported by iDutil Rovl (|l999l ). which for the less 
gas-rich simulations is —0.08 dex kpc~^ with a central oxy- 
gen abundance of 9.2. Thus, the initially simulated galax- 
ies a re located onto the observed mass-metallicity rela- 
tion (jTremonti et al.ll2QQ4l ). For the more gas-rich simula- 
tion (Sim VI), the initial metallicity gradient is somewhat 
steeper (—0.1 dex kpc~^) and the central metallicity lower 
(by 0.7 dex) than the less gas-rich galaxies. By doing that, 
these more gas-rich simulated galaxies are then initially onto 
the mass-meta ll icity relation observed at high redshift by 
iMaiolino et al.l (|2008l V Note that the exact adopted values 
are not critical since we are interested in the relative evolu- 
tion of the metallicity. 



3 METALLICITY EVOLUTION OF WET 
MERGERS IN THE LOCAL UNIVERSE 

In this section, we analyse the effects of galaxy interac- 
tions on the central metallicity in wet local encounters. 
Taking SimI as an example, we show in Fig. [1] the evo- 
lution of the gas-phase 0/H abundance of one member 
of the galaxy paii[j as the interaction proceeds. The evo- 
lution of the interaction can be followed by the relative 
distances between the center of masses of the interacting 
systems, plotted in the same figure. As it can be seen, 
there is a clear evidence of a decrease in the central gas- 
phase metallicity which can be correlated with the first and 
the second pericentre. Several previous works have shown 



The other member of the pair have a similar behaviour. 




100 ^ 



0.4 0.6 

time [Gyr] 



Figure 1. Evolution of the gas-phase oxygen abundance mea- 
sured inside a sphere of 2 kpc (black line), 4 kpc (blue line) and 6 
kpc (orange line) radius for one of the interacting galaxies of SimI. 
The relative distance between the center of mass of the simulated 
galaxies is also plotted (green line). 



that during a close encounter, strong tidal torques can de- 
velop, re-distributing angular momentu m and producing gas 
infall (e.^. [B arnes & Hernquist 199d;__ Mihos HernquistI 
'l996':'Tisseral l200q :,Di Matt eo et al.l l20Q8Vln the case that a 
metallicity gradient exists, the infalling gas will tend to have 
lower metallicity than that of baryons in the central region 
producing a dilution as it is clearly shown by our simulation 
(see also Rupke et al. 2010a; Montuori et al. 2010). 

In order to quantify the chemical changes as a function 
of time, we compute the central abundances inside concen- 
tric spheres of different radius in the range from 1 kpc to 
10 kpc. As an example. Fig. [1] shows the gas-phase oxygen 
abundance measured inside a sphere of 2 kpc, 4 kpc, and 
6 kpc of radius. As it can be appreciated, as one moves in- 
wards, the signals get noisier since the number of particle 
decreases. If one moves to larger radius, the signal decreases 
as chemical abundances start to be dominated by the exter- 
nal regions. As a compromised between these two effects, we 
select Rcen = 4 kpc (blue line) to define the central region 
where we will perform the analysis of central metallicities. 

In order to understand the decrease in the central abun- 
dances shown by Fig. [U we analyse the evolution of the 
central gas inflows and their abundances. We compute the 
central gas inflow at a certain time as the ratio between the 
new accreted gas mass within Rcen and the total gas con- 
tent within the same radius. Fig. [2^ shows the evolution of 
this central gas inflow as the interacting galaxies approach 
each other. We find that the gas inflow remains constant in 
around 2% until the first close passage when it rises up to 
^ 8%. From the second pericentre, the gas inflow steadily 
increases up to around 30% . This gas inflow is dominated 
by low- metallicity material as it can be appreciated from 
Fig. [2]3 (solid line) where we plot its 0/H abundance. The 
accretion of this low- metallicity gas dilutes the central oxy- 
gen abundance as shown by the continuous decrease detected 
after the first pericentre (Fig. [JJ), dotted line). The main 
drops in central abundances correlate with the occurrence 
of the pericentres driven by important gas inflows (Fig. [2^). 
Additionally, the gas inflow is also responsible for trigger- 
ing new star formation activity (Fig. [2]:). In SimI, the SFR 
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Figure 2. Evolution of gas inflow and of chemical abundances 
in the central region R < Rcen of one galaxy in SimI . a) Frac- 
tion of the gas falling within Rcen as a function of time. Dashed 
horizontal line represents the initial fraction, b) Oxygen abun- 
dance of the infalling gas (solid) and the mean central abundance 
within Rcen (dotted line), c) Evolution of star formation rate, d) 
a-enhancement of the infalling gas (solid) and its mean central 
value (dotted line). Arrows indicate the flrst and second pericen- 
tres. 



increases by a factor 3 from the first pericentre as the in- 
teraction continues fueling the gas inflow (Fig. [2^,c). As a 
consequence of this increase in the SFR, the enrichment by 
SNII also increases, diminishing the effects of chemical dilu- 
tion previously imprinted by the low-metallicity gas inflows. 

The increase of the SFR is also reflected in the a- 
enhancement, [0/Fe], of the gas in the central region (dotted 
line of Fig. [211). This central a-enhancement is present from 
the first pericentre but is much stronger around the second 
one. Note that the infalling gas into the central region is low 
a-enhanced as expected in a disc galaxy such as that of the 
Milky Way (solid line of Fig.[2]i), indicating that the increase 
of oxygen in the central region is principally due to SNII as- 
sociated to the central starburst. These general results ar e in 
agreement with those reported bv lMontuori et al ] ([2010') al- 
though these authors did not consider chemical enrichment 
self-consistently. As we will discuss in the next section, our 
simulations allow us to make f urther advances on th is issue . 

Previous works (e.^. Barnes HernquistI 1 19961 : 
iDi Matteo et al.ll2008l : ICox et al.ll2008l : iHopkins et al.ll200' " 



have already reported how star formation histories of 
interacting systems depend on the initial conditions and 
orbital parameters. Hence, the trends shown in Fig. [2] for 
SimI might vary among the analysed simulations as we will 
discuss in next sections. 



Figure 3. Evolution of central abundances as galaxy interaction 
proceeds for different simulations: I, II, IV and V. Red symbols 
identify the phase of interaction before the flrst pericentre, green 
points the phase between the flrst and second pericentres and 
blue ones the flnal merging phase after the second pericentre. 



3.1 Could the a-enhancement clock the 
interaction? 

As already shown for SimI, the a-enhancement starts to 
be detected after the first close passage. However, the most 
significant increase is recorded after the second one. In order 
to analyse if the g-enhanc ement could be used t o time the 
interaction as suggested bv lMontuori et aD (|2010l l, we study 
the ISM enrichment history along the interactions and its 
dependence on different orbital configurations. 

Fig. [3] shows the evolution of [0/Fe] in the central re- 
gions for different simulations with different orbital param- 
eters: SimI, Simll and SimV. We also include SimlV which 
assumes a lower energy per SN event but otherwise has the 
same parameters as Simll. We find that before the first close 
passage (red points) , the different simulations exhibit no sig- 
nificant variations in the 12 + log(0/H) nor [0/Fe] values. 
However, different feedback parameters and orbital configu- 
rations impact differently on the ISM chemical enrichment 
during the interacting phase between the first two pericen- 
tres (green points). As a general result, we find an increase of 
the [0/Fe], except for the coplanar retrograde orbit (SimV). 

In order to understand these differences in the a- 
enhancement. Fig. [4] shows the central a-enhancement and 
the SFR for one of the interacting galaxies in each simu- 
lation. As expected, the central a-enhancement (red dot- 
ted lines) correlates with the star formation activity (green 
lines) . We find that between the first and the second pericen- 
tre the [0/Fe] reaches the highest values for coplanar direct 
orbits (Simll and SimlV), decreases as inclination increases 
(SimI), and is almost negligible for the coplanar retrograde 
interact ion (SimV). Note that i n agreement with previous re- 
sults bv lDi Matteo et al.1 (|2007l l. we find that coplanar direct 
interactions trigger stronger starburts than coplanar retro- 
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grade configurations during the pre-merger phase (see Fig. 8 
of their paper). We also find that the simulation with stan- 
dard feedback (Simll) has slightly lower levels of the [0/Fe] 
and SFR, than the case of low energy feedback (SimlV) be- 
cause in the later, SN feedback drives weaker galactic out- 
flows, being less efficient at transporting material out of the 
central regions. 

As discussed previously by Montuori et al. (2010), the 
a-enhancement observed in interacting galaxies before the 
final coalescence of systems might date the interaction. How- 
ever, we find that this is valid when the a-enhancement is 
only regulated by the star formation process. In this case, 
the increase of [0/Fe] during the interacting phase would in- 
dicate the occurrence of a recent starburst (approximately 
< 1 Gyr), because the ISM have not had enough time to be 
heavily enriched by SNIa yet. On the contrary, the detec- 
tion of a decreasing [0/Fe] value would indicate that SNIa 
had time to release Fe into the ISM, which means that the 
time elapsed since the last starburst is longer than approx- 
imately ^ 1 Gyr. Our simulations which follow the dynam- 
ical and chemical evolution along the interaction suggest 
that other factors are also involved in the regulation of the 
a-enhancement. As can be seen for SimI, Simll and SimlV 
(Fig. 2]), the increase in the a-enhancement after the first 
pericentre is followed by a drop before the second pericen- 
tre. We find that this dilution of the central a-enhancement 
is caused by the low [0/Fe] enhancement of the very im- 
portant gas inflow (open blue symbols in Fig. |4]), triggered 
during the second close approach of the galaxies. This gas 
inflow vanish the clear traces in the central abundances left 
by the star formation activity, restarting the clock to date 
the interaction. 

Hence, our simulations suggest that the a-enhancement 
can only indicate the proximity to the first pericentre, pro- 
viding extra information on the parameters of the encounter 
since those with orbits close to coplanar retrograde configu- 
rations would tend to have very mild enhancement (SimV in 
Fig. 2]). However, our findings also suggest that the chance 
to observe close galaxy pairs with central a-enhancement 
might be low, because of the short time interval over which 
it could be detected (less than approximately 0.2 Gyr ac- 
cording to our experiments). 

3.2 Gas-phase metallicity profiles of wet 
low-redshift interactions 

As shown in the previous section, galaxy-galaxy interactions 
have a clear effect on the central abundances and the star 
formation of wet interacting systems at low redshift caused 
by the dynamical perturbation induced on the gas compo- 
nent. In order to illustrate how this affect the gas-phase 
metallicity profiles, we show a sequence of oxygen abundance 
maps which plot the spatial distribution for the interacting 
galaxies of SimI at different stage of the interaction, with 
their corresponding oxygen abundance profiles (Fig. [5|). It 
can be appreciated how the metallicity profile changes as 
the interacting system gets closer. This evolution is mainly 
related to the flattening of the profile caused by a decrease 
of central metallicities and a simultaneous increase of the 
abundances in the external region of the disc. We also show 
in Fig. [5] (right panel) that the abundance gradients com- 
puted by radial binning of all the gaseous particles are in 
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Figure 4. Evolution of a-enhancement of the central region (red 
dotted lines) and of the infalling gas (blue open symbols) and 
the star formation rate (green lines) for different simulations: I, 
II, IV and V. Arrows in each panel indicate the first and second 
pericentres. 



good agreements with the gradients determined with only 
the abundances of star-forming regions, as it is done with 
observations. 



3.2.1 Internal oxygen abundance profiles 

In order to quantify the evolution of internal metallicity pro- 
files, we carry out a linear regression fit to estimate the cen- 
tral oxygen abundances, (0/H)c, and slopes. The internal 
profiles are defined in the range of Sec — 7 kpc. The inner 
radial limit corresponds to three times the gravitational soft- 
ening used for baryons, and for the external radius, Re^t = 7 
kpc corresponds to approximately three scale-lengths of the 
initial stellar disc. 

Fig. [6] shows the linear fitting coefficients of the ra- 
dial profiles for different simulations. As it can be appre- 
ciated, there is a significant evolution of internal profiles 
after the first pericentre, which is consistent with a cen- 
tral dilution of abundances as expected from the discussion 
of the previous subsection. In agreem ent with observations 



OI tne previous subsection. In agreem ent witn observations 
(jKewlev et al.ll2006l: iRupke et al.|[2008l l and previous numer- 
ical works feu pke et aLl boiOaV the drop in (O/H)^ is about 
0.2 dex (red points). We also find that the final metallic- 
ity gradients are fiatter than those in earlier sta ges of the 
i ntera c tion (green p o ints) in agreement with iKewlev et al.l 
(|20i3);[RuEke et al.l (|2010l , and references therein). 

Note that although the evolution of internal profiles 
seems to be not very sensitive to different orbital config- 
urations in wet major mergers at low redshift, our results 
show some evidence that galaxies in coplanar direct orbits 
(Simll and SimlV of Fig. [6|) experience a larger flattening 
than those in coplanar retrograde ones (SimV in Fig. [6]). 
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Figure 5. Time sequence of the chemical properties in the gas component of Simll (before the first pericentre at t = 0.3Gyr (upper 
panel), after the first pericentre at t = O.SGyr (middle panel) and just before the second pericentre at t = 0.9Gyr (lower panel)) . Left 
panels: oxygen abundance maps. Middle panels: gas density map (in grey) and oxygen abundances of the star-forming regions (colour 
circles colour-coded as the left panels). Right panels: gas-phase oxygen abundance profiles. The black line shows the initial metallicity 
gradient. The magenta line and symbols show the gas-phase metallicity gradient computed with all the gaseous particles. For the top 
panels, only one galaxy is shown since the two galaxies are very distant and the field of view is smaller than the bottom panels. The 
right panels show that the abundance gradients computed by radial binning of all the gaseous particles are in good agreements with the 
gradients determined with only the abundances of star-forming regions, as it is done with observations. 



3.2.2 External oxygen abundance profiles 

As mentioned above, the external abundance profiles rise as 
the galaxy interaction proceeds. This is a common feature to 
all simulations, hence, we will take SimI again for illustration 
purposes. 

We define the external region of the gaseous disc as ev- 
erything outside i?ext = 7 kpc, which also denotes quite well 
the radius where the outer profiles depart from that defined 
by the central region (Fig. [5]). Then, we estimate the frac- 
tion of gas which has gotten into the external disc between 



two available time intervals, as well as the metallicity of this 
gas as a function of time. 

As can be seen from Fig. [Tt-, close to each pericentre, 
the external disc of interacting systems receives a significant 
fraction of gas, principally close to the second one (^ 30% 
at the second pericentre). This gas comes from the inner 
regions of the disc as the arms get opened and distorted by 
tidal forces (Fig. [5]). As a consequence, this gas tends to have 
high oxygen abundances (Fig.[7]3, solid line), which partially 
contribute, after the first pericentre, to increase the mean 
external abundances by ^ 0.4 dex (Fig. [7]3, dotted line). 
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Figure 6. Evolution of the linear regression coefficients ((O/H)^': 
red and Slope: green) of the oxygen abundance radial profile, 
computed in a range of Scq — 7 kpc for different simulations: I, 
II, IV and V. Arrows indicate the ffist and second pericentres. 

This result is in agreement with iRupke et all (|2010l ). who 
find that tidal tails transport metal- enriched gas from the 
galaxy inner regions to the outskirts. 

On the other hand, there is new star formation activity 
triggered in the external regions, on the spiral arms as can be 
seen from Fig. \5\ These new stars will also contribute with 
new SN II to the enrichment of the external ISM. Therefore, 
the evolution of the external abundances are produced by 
the combined effects of these two processes: i) the in situ 
SF activity and ii) the gas-rich material coming from inner 
regions as the interaction takes place and the spiral arms 
get distorted. 



4 METALLICITY EVOLUTION IN WET 
HIGH-REDSHIFT MERGERS 

In previous sections, we analyse wet low-redshift interac- 
tions of galaxies with around 20 per cent of the disc in the 
form of gas. Here, we will discuss the metallicity evolution 
in interacting systems when the initial amount of gas is in- 
creased up to 50 per cent (SimVI). This scenario is suitable 
for studying wet galaxy encounters at high redshift. 

For this simulation, we have chosen the same orbital 
configuration of Simll, but with a modified initial distribu- 
tion of baryons in order to start with more gas-rich systems 
(cf. Table [T|). We have also adopted the same star forma- 
tion and SN parameters as for Simll. The initial abundance 
gradient has been adjusted to be consistent with the obser- 
vations of the mass-metallicity relatio n at high redshift. Par - 
ticularly, we used the observations of iMaiolino et al.1 (|2008l ) 
at z ~ 3.5. Hence, the initial zero point is lower and the gra- 
dient steeper than in the case of less gas-rich interactions 
(see Section 2). 



Figure 7. a) Fraction of gas getting into the external region (R > 
Rext = 7 kpc) as a function of time, b) Oxygen abundance of the 
gas getting into this external region (solid line) and the mean 
external abundances (dotted line) as a function of time. Arrows 
indicate the first and second pericentres. These calculations are 
shown for one galaxy member of SimI as an example. 

Fig. [8^ shows the gas infiows for one of the galaxy mem- 
bers in SimVI. In this case, the amount of gas that falls 
into the central region is more significant than in its less 
gas-rich counterpart throughout the whole simulated time 
interval. Particularly after the second pericentre, the gas in- 
fiow roughly duplicates that of the less gas-rich interacting 
galaxies. 

The comparison of the SFR evolution of more gas-rich 
galaxies (Fig. [H})) with that of their less gas-rich counter- 
parts (Fig. [2]:) shows that long before the first close pas- 
sage, the SFR of the more gas-rich system is more active 
(around 30 M0/yr) than that of their less gas-rich counter- 
part. This is mainly due to the higher fraction of gas avail- 
able for st ar formation. Previous works (Bo urnaud et al.l 
l2007l . l2008l ) have suggested that gas-rich galaxies might ex- 
perience greater instability on the gaseous disc which can 
lead to an intense star formation activity via fragmentation 
and clump formation. In order to check this behaviour in 
our simulations, we compare the evolution of gaseous discs 
for galaxies in SimI and SimVI. Fig. [9] shows projected gas 
distributions of one member of the interacting system in 
SimI (upper panels) and in SimVI (lower panels), computed 
for different snapshots during the interaction. From these 
figures it is clear that, while the gaseous disc in the less gas- 
rich simulation evolves more smoothly, in the more gas-rich 
system it has highly fragmented and formed high density 
clumps where the gas satisfies the condition to form stars in 
agreement with previous workfl We checked that the Jeans 

^ Our less gas-rich simulation also shows clumps along the disc 
but with a significant lower level of fragmentation than in the 
more gas-rich counterpart. The contribution of these clumps to 
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mass is resolved by at least 4 gas particle so that fragmenta- 
tion can be numerically described. As it can be appreciated 
in the first low panel, the more gas-rich system shows high 
level of disk fragmentation well before the first close passage 
at 0.3 Gyr, when the two interacting galaxies are separated 
by more than 100 kpc and consistent with the peak in the 
SFR (see Fig. [SId). This suggests that the internal fragmen- 
tation of gas-rich discs is principally driven by their own 
self-gravity. To prove that this early strong SF activity is 
not induced by the interaction, we run an isolated galaxy 
case with the same star formation and SN parameters. The 
clump formation is also present in our isolated galaxy run. 
Previous works have shown that galaxy interactions rein- 
force the level of disc fragmentation for merg:ers at low and 
high redshifts (jTevssier et al.1 l2010l : iBournaud et al.1 l201lL 
respectively) . The detailed analysis of this issue will be part 
of a future paper focused on clump formation in isolated and 
interacting disc galaxies. 

After the intense SF activity reported during the early 
stage of interaction, the SFR in the more gas-rich galaxies 
begins to decline with two subsequent maxima shortly after 
the first and second pericentres, respectively (Fig. [SId). The 
induced star formation activity during the close encounters 
is less important than that expected considering the strong 
gas inflows experienced by these galaxies at these stages of 
evolution (Fig. [8^). We investigate which mechanism could 
be operating against the triggering action of gas inflows to 
inhibit the SF activity. We analysed the thermodynamical 
properties of the gas, finding that its internal energy (dotted 
line of Fig. [SId), increases continuously fed by the SN feed- 
back. This process is responsible of heating up and blowing 
away an important fraction of gas with an increasing effi- 
ciency as the galaxy approaches the second pericentre when 
it loses ^ 60% of its initial gas mass (^ 10% is lost at the 
first pericentre and more than 80% at the merging phase). 
As a consequence, the level of star formation activity is not 
only directly correlated with the strength of the gas inflows 
but it is also modulated by the SN feedback which is, in this 
simulation, more important than in the less gas-rich inter- 
actions because of the larger amount of available gas to feed 
star formation. 

We investigate the chemical evolution of the ISM by 
analysing the 0/H abundances and a-enhancements as in 
the less gas-rich interactions. Fig. [8]: shows the evolution 
of oxygen abundance of the infalling gas and of the central 
region. From this figure we can see that the central metal- 
licity of the gas increases steadily up to the first pericentre. 
We can also appreciate that the gas inflows, driven into the 
central regions during the first main starburst, are increas- 
ingly oxygen-rich as a consequence of the on-going impor- 
tant star formation activity. This is a consequence of the 
first early starburst which is not present in less gas-rich in- 
teractions. However, from the first to the second pericentre, 
the low metallicity of the infalling gas produces the oxy- 
gen dilution of the central region just as in wet local in- 
teractions (Fig. EJ)). The different levels of enrichment of 
these gas inflows are due to their different origins: high- 
metallicity inflows involved gas located nearby the central 



the SF activity in the less gas-rich mergers and a comparison with 
observations will be carried out in a future paper. 



region while lower metallicity inflows come from more exter- 
nal ones. Since, our discs have initially metallicity gradients, 
these different origins imply different levels of enrichment. 
To visualize this, in Fig. [10] (upper panel) we show the mini- 
mum, maximum, and the median location on the disc plane 
of those gas particles which will be immediately incorporated 
into the gas inflow^. Before the first pericentre most of the 
gas that falls into the central regions comes from about 4.5 
kpc and hence, is principally driven by internal dynamics 
while during the interactions the tidally driven gas inflows 
are able to transport material inward from further away re- 
gions of the disc as can be seen in Fig. [10] (upper panel). 

In Fig. [To] (lower panel) we show the stellar age gradi- 
ents as a function of time. We can see sharp drops in the pro- 
files as signatures of a discontinuous SF distribution along 
the disc, reflecting the clumpy SF activity. Profiles in blue, 
red, and green show the evolution of the SF distribution 
before the first pericentre. As shown by these profiles, the 
strong SF activity reported during this period (Fig. [Sj)) is 
located not only in the central region but also extended in 
the outskirt of the galaxy, contributing to enrich the ISM. 
As the interaction goes on, the main star formation activ- 
ity takes place principally in the central regions which gets 
rejuvenated. 

As a consequence of the contribution of a low- 
metallicity gas inflows during the close interaction (between 
the two first pericentres) the central 0/H abundance de- 
creases by ~ 0.4 dex (Fig. [8]:, dotted line). The fact that 
the central oxygen abundance decreases more slowly than 
that of the infalling gas, at this stage of the interaction, is 
a consequence of the SF activity being more concentrated 
in the central region, continuously supplying new chemi- 
cal elements via SNII events. After the second pericentre, 
we observe a second SFR enhancement driven by new gas 
inflow (Figs. [8t',b). This SF activity is very centralized in 
the nuclear region of the galaxy as shown in Fig. [10] (lower 
panel), producing the increase of the central 0/H abundance 
(Fig. [8]:, dotted line). 

Finally, we also investigate the a-enhancement of wet 
encounters at high redshift. Fig. [HJl shows that the initial 
[0/Fe] value is consistent with the a-enhancement observed 
in less gas-rich interactions (Fig.[2ll). The subsequent evolu- 
tion of the a-enhancement shows that immediately after the 
second pericentre and simultaneously to the SFR enhance- 
ment, the central [0/Fe] (Fig. [SJl, dotted line) exhibits a 
sharp increase produced by the important increase of SNII 
feedback associated to the more active central SF. After 
this a-enhancement, we find a significant drop in the central 
[0/Fe] values of around 0.8 dex. This drop is caused by the 
injection of Fe elements into the ISM by SNIa, whose pro- 
genitors formed ^ 1 Gyr before, during the main starbursts. 

4.1 Oxygen abundance profiles of wet interactions 
at high redshift 

In this section, we quantify the effect of wet interactions at 
high redshift on the gas-phase metallicity profiles. Similarly 

^ The location is determined as the distance of a given particle 
to the center of mass just before it moves into the central region, 
R<R cen • 
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Figure 9. Projected gas mass distribution onto the disc plane of one interacting galaxy for the less gas-rich simulation SimI (upper 
panels) and for its more gas-rich counterpart SimVI (bottom panels). Gas density maps are shown from the left to the right, at 0.3 Gyr 
when the galaxies are separated by ^100 kpc, well before the first pericentre; at 0.45 Gyr; just at the first pericentre at 0.7 Gyr and at 
0.8 Gyr between the first and second pericentres. 



to the case of wet interactions in the local Universe, the evo- 
lution of metallicity profiles in the more gas-rich interacting 
galaxies show a clear flattening and the central dilution of 
0/H abundances. We quantify the evolution of the metallic- 
ity profiles by performing linear regressions. Fig. [11] shows 
the (0/H)c and slopes as a function of time (red and green 
points, respectively). 

As expected, we find that the evolution of (0/H)c 
matches that of the central abundances computed within 
Rcen (see dotted line of Fig. \^). We can see an increase 
of the (0/H)c in around 0.6 dex during the approaching 
phase while after the first pericentre, the decrease reflects 
the impact of the low- metallicity inflows driven by the in- 
teractions. These low- metallicity inflows produce a decrease 
of ~ 0.5 dex in the oxygen central abundance while in less 
gas-rich counterparts, the decrease is of ~ 0.2 dex towards 
the flnal merging stage (Fig. [6]). 

The evolution of the gradients shows an approximately 
constant slope until the flrst pericentre where the proflles 
get more negative as a result of the oxygen carried into the 
central region by the early metal-rich inflows. However, from 
the flrst pericentre there is a continuous mean flattening 
of the metallicity proflles which gets even slightly positive 
around the second pericentre where again, it starts to be 
get steeper toward more negative values as more oxygen 
is pumped in by SNII produced during the second tidally- 
induced starburst (see Fig. [SId). 



5 DISCUSSION 

Several interesting aspects remain to be understood or even 
acknowledged in galaxy interactions and their impact on 



chemical properties. Numerical simulations are, indeed, the 
most suitable tool to study them. Different authors with 
different techniques and lev els of complexity have made im- 
portant contributions (e.g. IPerez etaP 1 20061 : iRupke et al.l 
l2010al : lMontuori et al.|[2Qlol ) to this area over the last years. 
All works agree to show the central abundance dilution due 
to low- metallicity inflows triggered during the interaction. 

Our simulations allow us to contribute further to ex- 
plore the role played by the induced SF and SN feedback in 
modulating the strength of this metallicity dilution and the 
evolution of the metallicity gradients. Particularly, we run 
Simlll which only follows the dynamics of the gas with flxed 
metallicity assigned according to an initial abundance pro- 
flle. No new star formation activity is permitted in Simlll. 
As it can be appreciated by comparing the evolution of the 
central oxygen abundances in Simll and Simlll in Fig. 1121 
neglecting the subsequent star formation and chemical en- 
richment leads to an overestimation of the metallicity dilu- 
tion (and flattening of the gradients). 

One of the effects of star formation is the generation 
of new SN II which contributes with new chemical elements 
and also pumps energy into the ISM. The strongest the gas 
inflows, the strongest the starburst s and consequently, the 
impact of SN feedback which can heat up and blow away 
signiflcant fractions of the remaining enriched gas, regulat- 
ing the subsequent star formation activity and modulating 
the chemical enrichment. In SimlV we use a lower energy 
input for SN event and thus, reducing the impact of energy 
feedback. As a consequence, a more intense star formation 
can develop producing a higher chemical enrichment than 
in Simll. From Fig. [12] we can see that the central oxygen 
abundance is very similar to that of Simll until the flrst 
pericentre while after that, the mean oxygen abundance in 
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Figure 8. High gas-rich simulation (SimVI): a) Fraction of gas 
inflowing within Rcen as a function of time, b) Evolution of star 
formation rate (solid line) and of the mean internal energy < u > 
of the gas (dotted line), c) Oxygen abundance of the infalling 
gas (solid line) and the mean central oxygen abundance (dotted 
line), computed within Rcen- d) a-enhancement of the infalling 
gas (solid line) and it mean central value (dotted line). Arrows 
indicate the location of flrst and second pericentres. 
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Figure 10. Upper panel: minimum (green), median (blue solid 
line) and maximum (red points) location on the disc of gas par- 
ticles that will be immediately incorporated into the gas inflow 
in the more gas-rich simulation (SimVI). Arrows indicate the lo- 
cation of flrst and second pericentres. Lower panel: stellar age 
proflles for different stages of evolution of SimVI, corresponding 
to: 0.3 Gyr (blue), 0.5 Gyr (red), 0.66 Gyr (green, flrst pericen- 
tre), 0.86 Gyr (magenta, after the flrst pericentre) and 1.16 Gyr 
(black, second pericentre). Vertical line indicates Rcen- 



SimlV is larger since the impact of the mass-loaded out- 
flows are weaker and are not so efficient at transporting ma- 
terial outside the central region and at inhibiting the star 
formation activity (the effects on the star formation rate, a- 
enhancement and metallicity profiles can be visualized com- 
paring results from Simll and SimlV in Figs. O [H and [6|). 
Hence, it is important to model the on- going star formation 
activity in a self- consistent way with energy and chemical SN 
feedback in order to have a more comprehensive picture of 
the evolution of the metallicity and gradients during galaxy 
interactions. 

As a consequence of this complex interplay of various 
physical processes, we find that some a priory obvious corre- 
lations are actually difficult to detect. One example of this is 
the use of a-enhancement to clock the interaction as we have 
already discussed in Section 3.1. Another expected correla- 
tion is that between the strength of the central starbursts 
and the metallicity dilution. In Fig. [T3]we plot the slopes of 
the metallicity gradients and the central oxygen abundances 
as a function of star formation rate for wet low-redshift in- 
teractions. We can see that globally there is a correlation so 
that the highest the star formation rate, the strongest the 




Figure 11. Evolution of the linear regression coeflicients 
((0/}i)c'- red and Slope: green) of the oxygen abundance pro- 
flles in the more gas-rich encounter (SimVI) computed in a range 
of Scg — 7 kpc. Arrows indicate the flrst and second pericentres. 
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time [Gyr] 

Figure 12. Effect of star formation and feedback: evolution of 
central gas-phase oxygen abundances (R < Rcen) of one galaxy 
in Simll (black full line), Simlll (orange dashed line), and Sim IV 
(red dash-dotted line). The relative distance between the centre 
of mass of the simulated galaxies is also plotted (green line). 

metallicity dilution and the shallower the abundance pro- 
files. However, the relations depend strongly on the orbital 
parameters. SimI and SimV have weaker starbursts than 
Simll but nevertheless, they experienced metallicity dilution 
in the central regions due to the tidally-induced gas inflows. 
Hence, according to our findings, it might be diflficult to es- 
tablish this correlation observationally where galaxy pairs 
can be observed at different stages of evolution and with a 
variety of orbital parameters. 

Regarding the gas-richness, from Fig llll we can see that, 
although there is a global trend in diluting central metallic- 
ity and in fiattening the mean gradients, there are also fiuc- 
tuations which correlate with the triggering of infiows and 
the injection of fresh oxygen by the new SN H. Nevertheless, 
the mean metallicity dilution is stronger than the less gas- 
rich counterpart (Fig.[6|). Hence, the gas-richness is another 
key factor. 

Depending on the astrophysical properties and orbital 
parameters, galaxies in pairs can be located differently on 
the mean mass-metallicity relation depending on the stage 
of the interaction at which they are observed. This effect 
might be particularly relevant for the study of high red- 
shift galaxies and could be at the origin of the apparent 
increase with redshift of the scatter in the mass-metallicity 
relation. It could also explain the presence of outliers in the 
local mass-metallicity relation. For example, more gas-rich 
systems observed just as they are first approaching might 
be dominated by the massively and clumpy star formation 
activity taking place at all radii, before the interaction can 
actually trigger larger scale infiows capable to transport low- 
met allicity gas into the central region. Hence, these galaxies 
migh t show an excess of oxygen (e.g. Michel- D ansae et al.. 
I2OO8I ). 



6 CONCLUSION 

By using hydrodynamical simulations of interacting pre- 
prepared galaxies of comparable masses, we study the chem- 
ical evolution of the interstellar medium during the interac- 
tions. Our simulations include self-consistently, SN feedback 



and a detailed chemical model, which allows us to follow the 
evolution of the gas-phase metallicity from the initial to the 
final stages of an interaction. These simulations constitute 
a unique set to study the effects of interactions on chemical 
properties. 

We analyse the different chemical response to galaxy 
interactions, varying orbital configurations, strength of en- 
ergy SN feedback and initial gas fractions in the galaxy discs. 
Particularly, we focused our study on wet interactions in the 
local Universe. For comparison, we also analysed a more gas- 
rich interaction in order to reproduce wet mergers typical of 
high redshift. 

Our main results can be summarized as follows: 

1) We find that a low- metallicity gas infiow developed 
from the first close passage dilutes the central oxygen abun- 
dance and triggers starbursts at the first and second peri- 
centre. We find that the amplitude of the central oxygen 
dilution of close galaxy pairs is fairly consistent with obser- 
vations. This dilution is observed with a large scatter for 
interacting systems separated less by than 30 kpc. The am- 
plitude varies according to orbital configurations and phases 
of interactions. 

During the interactions the metallicity gradients get 
quite shallower with fiuct nations, such as sometimes they are 
a slightly positive gradient. Both the release of new chemical 
elements and energy by SNII, produced during the tidally- 
induced starburst, modulates the strength of the central di- 
lution and the metallicity gradients. This behaviour seems 
in good agreeme nts with what it is observe d for galaxy pairs 
(observations bv lKewlev et al. I I2OO6I . I2OIQI ). We note that it 
is not possibl e from these simulations to test the claims of 
ICresci et al.1 (|201Q| 1 since our boundary conditions are not 
consistent with a cosmological scenario where gas accretion 
from the intergalactic medium can fuel the galaxies. 

2) We find an increase of the interstellar [0/Fe] val- 
ues, observed after the first pericentre as suggested by pre- 
vious works. This a-enhancement does not only refiect the 
tidally-induced star formation activity but is also affected by 
galactic outfiows and the abundances of the gas infiows. As 
a consequence, we show that the a-enhancement can only 
indicate the proximity to the first pericentre and in those 
encounters with non coplanar retrograde orbits. After the 
first pericentre, gas infiows might dilute the central [0/Fe] 
values, erasing the trace left by the star formation activity, 
and restarting the clock to date the interaction. We also no- 
tice that the period of time where this enhancement could 
be observed is quite short (< 0.2 Gyr). 

3) More gas-rich interactions show massive and clumpy 
star formation on the disc from early stages of the inter- 
action. During the interaction, from the first pericentre to 
the final stages, we find similar behaviours as in the case 
of the less gas-rich runs with the particularities imprinted 
by the more gas availability. Gas-richness has an impact on 
the strength of low- met allicity gas infiows which in the sim- 
ulated interaction, are able to dilute the central abundances 
by 0.5 dex, a factor of 2.5 larger than in its less gas-rich 
counterpart. There are important variations in the central 
abundances and gradients as the interaction of more gas-rich 
galaxies proceeds due to the competition between various 
processes, such as the injection of new chemical elements af- 
ter starbursts, the action of outfiows which blow away part 
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Figure 13. Metallicity gradients and central abundances as a function of star formation rate for interacting galaxies in SimI, Simll, 
SimlV, and SimV. 



of the enriched material, and the gas dynamics of clumpy 
discs. 

4) We find a correlations between the metallicity gra- 
dient, the central chemical abundance and the strength of 
the star formation activity. However, the slope and scat- 
ter of these correlations depend on orbital parameters, SN 
outflows and gas-richness. Hence, it is not a simple relation 
which can be easily unveiled observationally. The metallicity 
gradients are affected by the strength of the central nuclear 
starbursts and the extended star formation detected in the 
outer regions of the galaxies. Star formation in clumps are 
also found to play a role in modulating the metallicity gradi- 
ents. A detailed analysis of disc fragmentation as a function 
of gas richness will be carried out in a forthcoming paper. 

Our results suggest that metallicity evolution is closely 
related to gas dynamics during the encounters. Denying pro- 
cesses that affects the gas dynamics and its level of enrich- 
ment such as star formation and SN feedback, can contribute 
either to overestimate or underestimate the effects on the 
central abundance dilution or metallicity gradients. These 
findings show that there are still many issues to be addressed 
in galaxy interactions and their role in chemical evolution. 
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